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Abstract 

It has been suggested that a temperature gradient will induce a Leduc-Righi, or thermal Hall, 
current in the Majorana quasiparticles localized on the surface of class Dili topological insulators, 
and that the magnitude of this current can be related via an Einstein argument to a Hall-like energy 
flux induced by gravity. We critically examine this idea, and argue that the gravitational Hall 
effect is more complicated than its familiar analogue. A conventional Hall current is generated by a 
uniform electric field, but computing the flux from the gravitational Chern-Simons functional shows 
that gravitational field gradients — i.e. tidal forces — are needed to induce a energy-momentum 
flow. We relate the surface energy-momentum flux to a domain-wall gravitational anomaly via 
the Callan-Harvey inflow mechanism. We stress that the gauge invariance of the combined bulk- 
plus-boundary theory ensures that the current in the domain wall always experiences a "covariant" 
rather than "consistent" anomaly. We use this observation to confirm that the tidally induced 
energy-momentum current exactly accounts for the covariant gravitational anomaly in (1 + 1) 
dimensional domain-wall fermions. The same anomaly arises whether we write the Chern-Simons 
functional in terms of the Christoffel symbol or in terms of the the spin connection. 

PACS numbers: 73.43.-f, ll.10.Lm, 04.62+v, ll.15.Yc 
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I. INTRODUCTION 



One of the key properties of topological insulators is the intimate connection between 
the non-trivial bundle structure of the bulk electronic states and the presence of protected 
gapless surface modes. The most intuitive way of understanding this connection is that the 
twisted bundle gives rise to bulk quantum-Hall-like conductivities, and the gapless surface 
modes need to be present to soak up the corresponding conserved currents where they run 
into the surface of the sample Q, m this way the bulk-surface connection is seen to 
be a manifestation of the Callan-Harvey "anomaly inflow" mechanism [3(. Most of the 
Altland-Zirnbauer classes j^, 0] of topological insulators possess conserved U(l) charge or 
SU(2) spin currents, and the necessity of their protected surface modes can be understood 
via ordinary gauge-field anomalies. An important exception is the class Dili, which includes 
superconductors with spin-orbit interactions, and superfluid 3 He-B. Here the only conserved 
quantities are energy and (in the translation invariant superfluid) momentum. An anomaly- 
inflow understanding of the electrically neutral (Majorana) surface modes in the Dili systems 
therefore requires a failure of some edge-mode energy-momentum conservation law — in 
other words a gravitational anomaly jfj]. 

Gravitational anomalies originate in the A-genus contribution to the Dirac index theorem 
that is non-zero only in 4k space-time dimensions. They descend via a parity-violating grav- 
itational Chern-Simons term in 4k — 1 dimensions to an energy-momentum inflow anomaly 
in 4k — 2 space-time dimensions. For physically realizable topological insulators we are 
restricted to the k = 1, and therefore to a gravitational Chern-Simons term in a (2 + 1)- 



dimensional sur 



ace, and a gravitational anomaly in a (1 + l)-dimensional edge. 



Following |7Hll| we expect that (after the application of a small symmetry-breaking field 



that opens a gap) the (2 + l)-dimensional Majorana fermion surface modes of the Dili 
systems will include a gravitational Chern-Simons term in their low-energy effective action. 
It is argued in 12Hl4| that this term can, in principle, be observed through a thermal Hall (or 



Leduc-Righi) effect. A key step in the reasoning in [12Hl4j requires that, in analogy with the 
conventional Hall effect, a uniform gravitational field induces a surface energy-momentum 
current. The Leduc-Righi, coefficient is then obtained by means of an Einstein argument. 
The idea is that thermal equilibrium in the presence of a gravitational field requires the 
local temperature to vary so as to compensate for the gravitational red shift experienced 
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by radiation as it moves in the potential. The energy flux induced by a thermal gradient 
is then balanced by an equal and opposite energy flux due to the gravitational potential 
gradient. The thermal Hall conductance can thus be found from that of the gravitational 
Hall conductance. 



The purpose of this paper is to argue that, although the arguments in |l2Hl4J| are very 
appealing, the gravitational "Hall effect" is a little more complicated than its electromagnetic 
analogue. While a temperature gradient across a finite (2 + l)-dimensional surface does 
indeed induce a thermal Hall current whose magnitude is related to the gravitational anomaly 
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15], the surface-state energy gap exponentially suppress any surface thermal current. 



The heat must therefore be carried entirely by the (1 + l)-dimensional edge modes. In 
this respect the thermal current differs from the charge Hall effect, which can flow either 
at the (1 + l)-dimensional edge, or, in the presence of a uniform electric field, within the 
(2+l)-dimensional electron gas. Furthermore, the gravitational Chern-Simons term yields an 
energy-momentum flux that is proportional to to gradients of the Ricci tensor. Consequently 
a uniform bulk gravitational field cannot create an energy-momentum flux within the (2 + 1)- 
dimensional surface. A surface energy flux requires an inhomogeneous field — i.e. tidal 
forces. Nonetheless, the tide-induced energy-momentum flow does retain the bulk-boundary 
connection because it demands an anomalous (1 + l)-dimensional gapless mode to absorb 
the flux as it runs into an edge or domain wall. 

In section [TT] we will describe the thermal Hall effect and show how it can maintain an 
equilibrium balance between a temperature gradient and a gravitational potential gradient 
even in the absence of a bulk energy flow. In section II III we review the Callan-Harvey 
anomaly-inflow picture, and stress that this mechanism always leads to the covariant form 
of the associated anomaly. In section [IV] we explain the origin of the gravitational anomaly 
in (1 + 1) dimensional chiral theories. In section |V] we compute the energy momentum flows 
arising from a (2 + 1) -dimensional gravitational Chern-Simons functional and show that 
it exactly accounts for the anomaly obtained in section IIV1 We also show that the same 
anomaly is obtained from the Chern-Simons functional whther it is written in terms of the 
Christoffel symbols T or the spin connection u. Finally section IVT1 provides a brief summary 
of our results. 
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II. THERMAL HALL CURRENTS 



The edge of a (2+l)-dimensional quantum Hall system hosts gapless chiral fermions 



161 ] , and both the edge of a p x + ip y superconductor and domain walls on the surface of a 



suitably engineered topological insulator host chiral Majorana fermions [ill, |17|]. Suppose 
for a moment that these modes can be modelled as a set of n independent conformal fields 
having (positive or negative) propagation velocities i — 1, . . . , n and central charges q. 
Then, at temperature T, each independent chiral edge-mode contains an energy density [IS] 

~ - ~ 77 ^-T 2 , (1) 



12 k h 



where ks is the Boltzmann constant. Thermal energy is therefore being transported along 
the edge at a rate [3] 



ViSi 



= ^( c _c)^T 2 . (2) 
12 v 7 v 7 

Here c and c are the total conformal charges of the right and left-moving modes respectively. 
Although motivated by the model of independent modes, this formula continues to hold for 
more complicated conformal edge-mode theories [15 ] . If we construct a parallel-sided Hall 
bar and maintain a small temperature gradient AT across it, then the difference between 
the contra-propagating energy fluxes (|2J) on the two edges gives rise to a net thermal current 

Jl-r = Cl-rAT (3) 

that flows along the bar and perpendicular to the temperature gradient. Here 

C L _ R =(c-c)^T, (4) 

is the Leduc-Righi coefficient. 

It is remarkable that the non-universal edge-mode velocities have cancelled, leaving in 
Cl-r/T only fundamental constants and the numbers c, c that are characteristic of the 
quantum Hall phase. It is therefore reasonable to suppose that Cl-r/T may be extracted 
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from topological data, as is the quantum Hall coefficient. It is, however, difficult to provide 
a direct derivation of thermal conductivities from linear response theory. There is no term 
that can be added to the Hamiltonian to describe the temperature. An ingenious trick 
was introduced by Luttinger [^(J, who instead coupled the system to gravity and proceeded 
indirectly by adopting the method used by Einstein to relate diffusion coefficients to viscosity 
2l| . Luttinger's idea is that the de-equilibrating effect of a small temperature gradient will 
be precisely compensated for by a gravitational potential $ provided that 

IdT 1 <9$ 



T dx cjL t dx 



(5) 



Consequently, assuming that all currents vanish in equilibrium, and that the effects of the 
two driving forces are additive, a linear-response derivation of the current induced by gravity 
allows one to find the current induced by the thermal gradient. 

In the present case imagine a rectangular Hall bar whose upper, right-propagating, edge 
at co-ordinate y — y\ is held at temperature 7\ and whose left-propagating lower edge at 
y = yi is held at temperature T 2 > 7\. If the bar lies in a gravitational field such that the 
gravitational frequency shift obeys 



v(yi) _ / £00(3/2) = Tx 
u(v%) ~ V S'oo(yi) T 2 ' 

then, as the thermal excitations from the hotter lower edge rise on the left hand vertical side 

to the upper edge they will red-shift to the lower temperature. Similarly, as the excitations 

from the cooler upper edge descend via the right hand vertical side to the lower edge they 

will blue-shift to match the hotter temperature. The system is in equilibrium therefore. 

Since for weak gravitational fields we have 

vS^ + T^, (7) 

C light 

this situation satisfies ([5]). Observe, however, that in our Hall bar, the currents are not zero 
in equilibrium. Therefore a knowledge of the thermal Hall current does not allow one to 
deduce the gravitational Hall current, nor vice versa. 

The steady-state equilibrium of the Hall bar requires no thermal or other form of energy 
to be flowing within the gapped surface states. Indeed an extra energy-momentum flux into 
the gapless states on the vertical sides would mess things up. This already suggests that a 
uniform gravitational field should not induce a surface energy flow. 



This suggestion is perhaps not surprising. A mathematical analogy between the con- 
ventional Hall effect and gravitation would naturally identify the field strength F with the 
Riemann curvature R. A uniform field gravitational field does not, however, require space- 
time curvature. The Rindler metric 

dr 2 = ( 1 + ^plA d t 2 - J-dr 2 , g = cf ight /r (8) 

\ C light J C light 

of a uniformly accelerated observer provides a gravitational potential $(r) = (r — r )g, but 
is merely a re-parametrization 

r i / c light * 

CiightT = r smh ' 
X = r cosh 
of a part of Minkowski space with flat metric 



r 

( Clight M 
V r o ) 
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dr 2 = dT 2 - -^—dX 2 . (9) 

C light 



It might, therefore, be more physical to identify the thermal-Hall analogue of the electric 
field with the Christoffel symbols V which describe the frame- dependent inertial forces that 
we perceive as gravity. If this new analogy is to work, the energy-momentum influx into 
the edge modes would have to be given by the non-covariant "consistent" gravitational 
anomaly, which contains T's, rather that the "covariant" anomaly which contains only R 
22[. In the following sections, however, we will argue that anomaly inflow always give rise 
to the covariant anomaly, and not to the consistent anomaly. Moreover, we will see that 
gradients of curvature are needed to produce an energy flow into edge states. 



To simplify the argument, we will follow the authors of [12| and argue that since we are 
interested in topological effects, we can choose non-universal quantities such as propagation 
velocities as we like. We will therefore from now on make all modes propagate at the speed 
of light, work with fully relativistic systems, and use natural units, in which h = cu g ht = 1- 



III. THE CALLAN-HARVEY MECHANISM AND COVARIANT VERSUS CON- 
SISTENT ANOMALIES 

Let us recall how the conservation (or non-conservation) of a gauge current is related 
to the gauge invariance (or the lack of it) of an action functional. Suppose, for example, 
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that S[A] is a functional of an su(N) Lie-algebra-valued gauge field = \ a A®, where 
the matrices A a are the generators of su(N). We define the matrix- valued gauge current 
J»(x) = \ a J^ a by setting 

5S[A} = J d d xtr{J^5A^}. (10) 

Under a gauge transformation the field changes as A^ — > A 9 ^ = g^A^g + g^d^g, where 
g G SU(iV). For an infinitesimal transformation g — 1 — e the transformation becomes 
A^ — > Ap + SeAp, where <5 e A M = —([A^, e] +d^e) = — V M e. The corresponding change in S[A] 



is 



6 e S = J d d xtr{J»([e,AJ-d^)} 

= J d d xtr{e(d^ + [A^,J^)}. (11) 



The covariant divergence 



V^ = d^ + [A„J^} (12) 

is therefore zero if and only if S[A] is gauge invariant. 

We are interested in effective actions S[A] that arise as a result integrating out a collection 
of Fermi fields ip, ^ in the presence of a classical background gauge field A^. 

eM~S[A]} = J d[#f[^]exp{-%;,VM} . (13) 

The calculated currents are then the expectation value J M = (J^) of a quantum operator. 
The original S[ip, ip^ , A] action will be invariant under A^ — > A 9 ^, ip — > g^ip, ^ — > ip^g, but 
the invariance may lost during the functional integration. In this case we will have 

= G(A), (14) 

where the anomaly G(A) is a local polynomial in the A^ and their derivatives. A gauge 
anomaly provides an obstruction to a subsequent quantization of the A^ fields, but when 
the A^ are simply classical probes it provides a useful source of non-perturbative information. 

The Callan-Harvey effect links the non-conservation of gauge and other currents to an 
inflow of charge from some higher dimensional space in which the anomalous theory is 
embedded as modes localized on a domain wall or string defect. In the cases we are interested 
in, the inflow is derived from a Chern-Simons term in one-higher space dimension. 
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As usual we will think of A as a Lie-algebra- valued one- form A = \ a A a ^dx^, and define 
the field strength as the Lie-algebra-valued two-form 

F = dA + A 2 = -F^dx^dx" . (15) 

The Chern-Simons form u% n -x{A) is then defined as 

u 2n ^(A) = n I tviAF?- 1 }, (16) 
Jo 

where F t = tF + t(t — 1)A 2 . It is constructed so that du 2n -i = tr {F n }. For example, 

u 3 (A) = tr {AdA+^A 3 }, 

= tr {AF-^A 3 }, (17) 



and 



u 5 (A) = tr {A(dAf + ^A 3 dA + ^A 5 } 

= tr {AF 2 - X -FA 3 + ^A 5 }. (18) 



The F-free last term oc A in the second forms of to 2n -i has coefficient 

r -( ir - i n!(n ~ 1)! 

n_l J (2n-l)!" 



(19) 



It is this last term that governs the change in integrals of u 2n -i under large gauge transfor- 
mations. If A undergoes a finite gauge transformation 

A ^ A 9 = g~ x Ag + g~ x dg, (20) 

then 

U2n-i(A 9 ) = u^A) + c n ti {(g^dg) 2 ^ 1 } + da 2n _ 2 (A,g), (21) 
where, for example 25] 

a 2 = -tr{dgg~ l A} (22) 

and 

a 4 (A,g) = ~tr{{dgg- l ){AdA + dAA + A 3 ) - l -{dgg- l )A{dgg- l )A - {dgg- l ) 3 A}. (23) 
The Chern-Simons functional C[A] is defined by setting 

c ^^{^)^J u ^- Ml (24) 



where M is some 2n — 1 dimensional manifold. The coefficient in front of the integral has 
been chosen so that exp{iC[A]} is single- valued when M is the (2n — l)-sphere. In this case 

CW - C[A] = 2, (±J |-2| / tr (25) 



and it is shown in [24| that the right-hand side of (125]) is 27r times an integer whenever 
g G GL(n, C) or any of its compact subgroups such as SU(iV). This means that when a 
Chern-Simons functional appears in a functional integral 

Z = J d[A] exp{ikC[A] + ■■■} (26) 

then gauge invariance demands that k be an integer. This constraint on k need not hold 
when C[A] appears in an effective action. Indeed k is 1/2 when we integrate out a massive 
Dirac fermion in odd dimensional space time. 

Given a (2n — l)-manifold M possessing a 2n - 2-dimensional boundary dM, we can 
use C[A] to construct an action S[A, g] = C[A g ] that is obviously invariant under A — > A , 
g — > hr x g. In this action, the gauge non-invariance of the bulk Chern-Simons term C[A] is 
compensated by the complementary gauge non-invariance of the Wess-Zumino action 

W[A,g] % C[A*\-C[A] 

= 2tt (j-J 1 [J a2n - 2 (A,g) + c n £tr{(rt) 2 "- 1 }} . (27) 

Although requires g to be defined on the 2n — l dimensional manifold M, the identity 

(MGr 1 ^)* 1 " 1 } = (2n - l)dtr {(p- 1 ^)^" 1 ^) 2 "- 2 } (28) 

ensures that variation of VK[A, gr] depends only on the values that 8 A and 5g take on the 
boundary dM. It can therefore serve as an anomaly- capturing non-local effective action for 



a 2n — 2 dimensional theory [23|. The meaning of the gauge-group element g depends on 
the context. In a two dimensional boundary g(x,t) could be the dynamical chiral boson 
equivalent to a chiral fermion. In two or higher dimensions it might parametrize a Higgs 
field that gives a left-handed chiral fermion a mass by coupling it to a right handed chiral 
fermion that does not itself couple to A. 

The gauge anomaly in a Wess-Zumino action for a four dimensional theory may be read 
off from 



/ d 4 xtr{eV,J^ z } = S t W[A,g] 

JdM 



S e C[A] 

l — [ tr {de(AdA + dAA + A 3 )} 
wr J aM 



24vr 
1 

24^ 



^ / tr {ed^A^Ar + d v A a A T + A^ CT A r )}e^ CTT d 4 x.(29) 



So 



tr {eV M J^ z } = -Ltr {ed^A„d a A T + d v A a A T + A^A^)}^. (30) 



Because this anomaly is found as the variation of the functional VK[A g], it satisfies the 
Wess-Zumino consistency condition 

It is therefore known as a "consistent" anomaly. The right hand side of the (non) conserva- 
tion equation is not gauge covariant, however, and so neither is the left. The gauge current 
itself is therefore not covariant, and the physical meaning of the (non) conservation equation 
is unclear. 

In the full bulk-plus-boundary theory, whose gauge-invariant effective action is C[A 9 ] the 
non- zero divergence of the boundary current is being supplied by the inflow of gauge current 
from the higher dimensional bulk. This bulk current is covariant, 

tr {X a J x } = -Ltr {\ a F, v F aT }e x ^ T . (31) 

OZ7T z 



It comes from the variation 
8 / u 5 



3 I tr{5AF 2 } + [ tr {5A(AdA + dAA + -A 3 )} 

JM JdM 2 

= 3 I tr{5AF 2 } + [ ti{5A{AF + FA-\a 3 )}. (32) 

JM JdM 2 

We usually ignore the boundary term when computing a bulk current, but in the total bulk- 
plus-boundary theory we must retain it as it provides a contribution to the current in the 
boundary of 

tr {\ a X»} = -Ltr {X a (A u F aT + F ua A T - A v A a A T )}e^ . (33) 

This quantity is exactly the extra current ([^ equation (2.16)) that has to be added to the 
consistent current to obtain the covariant anomaly 

tr {A a V M (J^ z + X»)} = ^tr {\ a F^F UT }e^ T . (34) 
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The new current Jfr ot = J^ z + is now gauge-covariant, and its anomalous divergence 



entirely accounted for by the Callan-Havey anomaly inflow [26| 
Similarly, in two dimensions we find that 



IB- 



is the consistent anomaly, and 



V,^ z = ^»d,A„ (35) 



X m = 1-^A V (36) 

4:71 



is the Chern-Simons term's contribution to the boundary current. Then 
V^J&z + X") = ^d,A u + ^{d,A v +[A^A v }) 

47T 47T 

= ^{d^-d^ + iA^AA) 

= ^f F *» (37) 

is the covariant anomaly. 

We have seen that Bardeen-Zumino polynomial X^(A) that converts the consistent gauge 
current to the covariant gauge current is precisely the contribution to the boundary current 
provided by the boundary variation of the bulk Chern-Simons functional. The analogous 
conversion of a consistent to a covariant gravitational anomaly requires an extra integration 
by parts, and so is more intricate. Indeed some puzzlement was expressed in js| about what 
happened to the inflowing energy- momentum — see the discussion after equation (30) in {3] 
— but it was later understood that the anomaly inflow is always leads to a covariant current 

In the above examples, the Chern-Simons term was defined in the bulk and the lower- 
dimensional degrees of freedom resided on the boundary. This is, for example, the situation 
in the ordinary quantum Hall effect. For (3 + l)-dimensional topological insulators it is 
the Chern-Simons functional that is defined on the boundary, and the lower-dimensional 
theory is defined on a domain wall within the boundary. In this case the coefficient of the 
Chern-Simons functional is multiplied by sgn(m)/2, where m denotes the mass gap induced 
by a small symmetry breaking field that changes sign at the domain wall. The resulting 
domain-wall chiral fermions then experience half of the the usual inflow from each side, but 
there are two sides, and so the resulting edge-theory anomaly is unchanged. 
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IV. TWO-DIMENSIONAL GRAVITATIONAL ANOMALIES 



In this section we will review the origin and possible forms of gravitational anomalies. We 
start from an effective action S[g] that depends on the space-time metric g^ v . The associated 
Hilbert energy-momentum tensor T^ v is then defined by the variation 

SS eS = - l -J d d x^T^5g, V) (38) 
= +\f d d x^T, u 5g^. (39) 



2 . 

Under a change of co-ordinates x M — > x'^ = x^ + e M we have g^ u — > g' = g^ v + Sg^ u , where 

5g^ = {C t g)^ 

= + V„e„. (40) 

Here C e g denotes the Lie derivative of the metric with respect to e M , and V M is the covariant 
derivative with respect to the torsion-free Levi-Civita connection. When the effective action 
is invariant under this reparametrization we find (taking into account that J 7 *" = T VfM ) that 

= - J d d x^T^V,e u 

= J d d x^e u V,T^. (41) 

Thus a gravitational anomaly — i.e. a failure of the covariant conservation law V M T^ = - 
reflects a failure of reparametrization invariance. While it seems reasonable that any physical 
system should be independent of how we choose to describe it, co-ordinate dependence can 
creep into S[g] when we tacitly tie a regularization procedure to the co-ordinate grid rather 
than to some intrinsic property such as the metric. 

An equivalent Lorentz anomaly can also occur in theories when we use a frame field 
rather than the metric to encode the geometry. This anomaly manifests itself as a failure of 
the energy momentum tensor (now defined in terms of a functional derivative with respect 
to e£) to be symmetric. 

We will focus on two-dimensional systems expressed in terms of Euclidean signature 
isothermal co-ordinates x, y, in which ds 2 = e^(dx 2 + dy 2 ). It is convenient to set z = x + iy, 
z = x — iy so that ds 2 = e^dzdz. The non-zero component of the metric tensor and its 
inverse are then g Sz = g z2 = (l/2)e^, and g zz = g zz = 2e~*. In these complex isothermal 
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co-ordinates the only non-zero entries in the Levi-Civita connection are 

I* = d- z <f>. (42) 

The curvature is completely encoded in the Ricci scalar 

R = W\ v = R zz - zz + R z \- Z = -4e"*C-0- (43) 

In our convention R is twice the Gaussian curvature, and hence positive for a sphere. 

The effective action for a left-right symmetric theory with conformal central charge c was 
obtained by Polyakov 



28( as 

SpolyakovM = / d 2 x(8(j)f 

yo7r J 



c 



o} _ , d 2 xd z <pd- z <p. (44) 

Here d 2 x denotes dx A dy = dz A dzjli. To evaluate (|4"4"I) for a given geometry we must 
select a system of isothermal co-ordinates, and this choice is not unique. It is therefore not 
immediately obvious that SWyakov [#] is co-ordinate independent. To verify that it is so, we 
must examine the conservation of the energy-momentum tensor. 

Now to make use of the Hilbert definition of T^ v , we must be free to make an arbitrary 
infinitesimal variation in the metric. A general variation, however, will take us away from 
the class of isothermal metrics. We therefore make a variation Sg^ u and follow if with a 
change of co-ordinates 

z — > z' = z + e(z, z) 

z -> z' — z + e(z, z) (45) 

so as to return to the isothermal gauge. Now 

5(ds 2 ) = [e' l '(ed z (f) + ed z (f) + d z e + d z e) + 5g zz + 5g z2 )]dzdz 

+(Sg zz + e^d z e)dzdz + (5g zz + e*d g e)dzdz. (46) 

The required co-ordinate change is obtained by solving 

e+d,e = -6g„ 

e<f>d- z e = -5g- zz . (47) 
13 



Let us assume for the moment that given 5g zz and 5g zz we can always solve these equations 
for e and e. Then, comparing with 5(ds 2 ) = e^5<f)dzdz we find that the metric variation leads 
to 

6<j) = ed z (f) + ed-^ + 8 z e + d- z e + e~*{8g- zz + 5g zz ) . (48) 

We insert this variation of 4> into equation (144p . and, assuming that integration by parts is 
legitimate, reduce the terms involving e to 



c 



127T 

c 
12tt 



J d 2 xd- z e Q(^0) 2 -<9 2 2 >) 

J d 2 xe~Hg zz (dlJ - \(d z <P)^j . (49) 



On comparing with 



^SpoiyakovM = ~\ J d 2 x ^/gSg^T^ 

d 2 x^5g- z - z T ss , (50) 



where ^fgd 2 x = e^dxdy, we read off that 



Thus 



Similarly we find that 



9 ZZ 9 ZZ T ZZ 

Ae~ 2 *T zz . (51) 



T - = ^{ d ^-\^) 2 )- (52) 



Next, examining the effects of + 5g zz , we have 

^ Polyakov [(7] = j d 2 xe-^6g- zz + 5g zz )(-d 2 - z <f>). (54) 



From this we read off that 



6n 



' ; e" 2 ^ 2 0, (55) 



247T 
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and 



We also recover the well-known trace anomaly 29 j 



T% = g- zz T~ zz + 9z-zT" z = e*T zS = ^-R. (57) 
This is a comforting consistency check, as Polyakov derived f|44|) by working backwards from 

We can now verify that T^ v is covariantly conserved: 

^(vr„ + v%.) = v,r„ + v,r B 

w z ± zz ~ w z ± zz 1 zz- 1 zz 

= d- z T zz + d z T- zz - d z 0T- zz 

= 0. (58) 

This is evidence that Sp i y akov [<?] is indeed co-ordinate independent. There is a problem how- 
ever: if Spoiyakov is co-ordinate independent then its functional derivative T^ v must transform 
as a tensor. When we make a holomorphic change of variables z = z(£), z = z(Q, however, 
we have ds 2 = e x d(d( = e^dzdz and so 

Consequently 



+k i'- (59) 



24n \dz J V 2 J 24tt \ C' 2 \C 

where is the energy- momentum tensor component evaluated in the (, ( co-ordinates and 
{(, z} is the Schwarzian derivative in whose definition the primes denote differentiation with 
respect to z. Our T^ u does not transform as a tensor therefore. The paradox is resolved by 
looking back at the first line in equation ( H9|) . We see that if we are allowed to integrate by 
parts we can take the d z derivative off of e and onto T zz . Thus any holomorphic addition 
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to T zz is invisible to the variation 5g zz . Another way of saying this is that there can be 
metric variations 5g zz that cannot be written in the form the form Sg zz = —d z e = — 2V«ea. 
(The displacements e and e should really be written as e z and e z as they are the components 
of a contravariant vector.) The solvability of (|47l) depends on the global topology or on 
boundary conditions. On a torus, for example, metric variations due to change in the 
modular parameter r are not expressible in this way. On a closed manifold of genus g > 2, 
there will be 3(g — 1) linearly independent unobtainable metric variations. 

The addition of a purely holomorphic term is indeed required. The full operator energy 
momentum tensor is 

tz = T(z) + ^-^ 2 0-i(^0) 2 



C ^2 



24vr 



(61) 



where, for a free c = 1 boson field <p(z, z) = <p(z) + <p(z) for example, 

f(z) = : d t <p(z)dzV>(z) : 

= lim (d z ip(z + 8/2)d t ip{* - S/2) + . (62) 

(Note that conformal field theory papers often define T(z) to be — 2it times fl62l) so as to 
simplify the operator product expansion.) The operator T(z) has been constructed to be 
explicitly holomorphic, but at a price of tying its definition to the z, z coordinate system - 
both in the mode normal ordering expression in the first line and by the explicit counterterm 
in the second. It is not surprising, therefore, that under a holomorphic change of co-ordinates 
the operator T(z) does not transform as a tensor. It is well-known that instead 

f W = (g) 2 f(C)-^{C,*}. (63) 

We see that the inhomogeneous Schwarzian derivative terms cancel in the transformation of 
the energy momentum tensor defined in (161I) . Thus T^ v transforms as a tensor and is 
still covariantly conserved. It is notable that both the covariant conservation and the trace 
anomaly in are accounted for by the c-number terms. These properties are therefore 
independent of the quantum state in which the expectation is taken. This quantum state 
only influences the holomorphic part of (T zz ) and an antiholomorphic part of (T zz ). 
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In a chiral theory we might constrain both T 2z and T zz to be zero, while keeping the 
covariant form of T zz defined in the first line of (161"]) . The term V Z T ZZ needed for the 
continued mathematical validity of (158]) would then be interpreted as 

' d z R, (64) 



48vr 



so that conservation law (1551) is reinterpreted as the anomaly equation appearing in 15 1 

V z f zz = ~-^-d z R. (65) 

487T 

By adding in an identically zero term we can write this as 

V z f zz + V~ z f- zz = --^-d z R, (66) 

407T 

which at first glance looks like a covariant tensor equation. It is is not, however, because 
replacing the free index z with z leads to 

V z f z - Z + V~ Z T- Z - Z = --^-d- z R (67) 

on which the left hand side is identically zero, but the right need not be. Thus (166|) is not 
the covariant anomaly. 



A more symmetric treatment 33J divides the trace anomaly between the left and right 



chiral sectors and constrains one of them to zero. Then T zz remains zero, but 



so that 



f» u = —R. (69) 
M 48vr V ; 



This physical reinterpretation makes the (still mathematically valid) equation (155]) read 

V z f zz + V~ z f- zz = ~^-d z R, (70) 

V z f zz + V~ z f zz = +-^-d g R, (71) 
yo7r 

where the second term on the left hand side of the second equation is constrained to be zero. 
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In our z, z co-ordinates system we have y/g = \J—gzz9zz = —ie^/2, and we can write 
these last two equations in a covariant manner as 

V z fzz + V~ z f- zz = i-^-^ez-zcfR, (72) 
V z f z - Z + V~ z f zz = t^^ge- zz d z R. (73) 

907T 



In general euclidean co-ordinates we therefore have 



24 1 



= i-^y/ge va dPR. (74) 

The factor "i" appears in (J75l) because it is only the imaginary part of the Euclidean effective 
action that can be anomalous jf], [i(3] . It is absent when we write the equation in Minkowksi 
signature space-time, where 

V M T^ = JLJ-f'drR. (75) 
96-K^/g 

Note that f J75|) can be rewritten as V^T^ = where 

fnv = fp, _ _^ — L e ™R. (76) 

The new tensor T^ u is conserved, but not symmetric. We have therefore exchanged a 
reparametrization anomaly for a Lorentz anomaly. 

We now show that the manifestly covariant anomaly ( I75|) is that expected from the 
anomaly inflow. 



V. GRAVITATIONAL CHERN-SIMONS TERMS 

In this section we will use both the co-ordinate and frame-field (vielbein) description of 
geometric quantities. Thus e£ are the components of the frame field e a = e^d^ and e* b the 
components of the co- frame e* a = e^dx^, with 5% = e* a e^. The frame metric 

Vab = g^e&l (77) 

is diag(l, 1, 1) and diag(l, —1, —1) in Euclidean and Minkowski space, respectively. 

A gravitational (2 + 1) dimensional Chern-Simons functional can be written either in 
terms of the Christoffel-symbol form Y^ v = T^ ua dx a as 

c[T} = ^- [ tr{rrfr + |r 3 }. (78) 

967r J M 3 
18 



or in terms of the spin connection u a b = ^b^dx^ as 

CM= ^I tr{w<iw+ r 3} - (79) 

The integrands in these two functionals have the same exterior derivative 

dtT{TdT + ^T 3 } = dti{uduj + ^co 3 } =tr{R 2 }, (80) 

and so they coincide when M = dN is a boundary, but they are no longer equal when M 
itself has a boundary. Their normalization is related to the index 

Index(D Dirac ) = DimKer(D Dirac ) - DimKer(Dj )irac ) 

= i^X tr{fl2} (81) 

of the four dimensional Dirac operator. The Dirac index an even integer for any four- 
dimensional manifold possessing a spin structure. 

The spin connection is related to the Christoffel form by a GL(3) gauge transformation 

"V = <TV£ + efdtf, (82) 

and so 

C[u>] = C[T] - -1- / tr{(rfee*)r} - / tr{(e*rfe) 3 }. (83) 

96vr J aM 2887T J M 

Here the matrix-valued one-forms dee* and e*de are defined by (dee*Y u = (d a e%)el a dx a and 
(e*cfe) a 6 = e^d a e^dx a . 

The functional C[r] is invariant under reparametrization x M — > X^(x) up to boundary 
terms. To obtain an energy-momentum conserving theory it has to be attached to a suit- 
able boundary theory with compensating transformation properties. We do not have to 
write down the corresponding Wess-Zumino action W(T,X) to know the boundary theory 
anomaly. All we need to do is calculate the out-flowing bulk energy-momentum flux by 
computing the response of C[T] to a change in the metric. 

The variation of the Chern-Simons functional due to a change in T is 

sclT]= ^L ti{sm+ ^L ti{STir] - (84) 

To compute the contribution to the energy-momentum tensor we also need 

5r" w = ^{VJgxa + V a 8g aX - V x 5g V a). (85) 
19 



Then, making use of properties of the Riemann tensor that are unique to three dimensions 



(See |3l|, |32] for more details), we find 

SC[T) = -|- / Sx-sfgC^bg^ + boundary terms, (86) 

4o7T J M 



where 

= {e^VpK + e^V p <) (87) 

9 

is the Cotton tensor. We read-off the bulk energy-momentum tensor to be 

T „u = £_ 

24tt v ; 

In deriving this result we have had to integrate by parts a second time so as to remove the 
derivatives from the metric variations. Consequently the boundary terms are more compli- 
cated than the usual ones arising from the variation of gauge field Chern-Simons functionals. 
We are, however, confident that these boundary terms provide the same conversion of the 
consistent anomaly of the boundary theory into the covariant anomaly that we saw with the 
gauge anomalies. 

We restrict ourselves to product metrics of the form 

ds 2 = (dx 2 ) 2 + g ab (x°, x 1 )dx a dx\ a, b = 0, 1 (89) 

with boundary being at x 2 = 0. The Ricci tensor apearing in ( IHTl) then coincides with the 
Ricci tensor of the two-dimensional boundary, and can be written as 

R% = ^R(x°,x 1 ), a, b = 0,1 (90) 

The flux of the a = 0, 1 energy-momentum components into the boundary becomes 

The energy momentum inflow into the boundary therefore precisely accounts for the the 
gravitational anomaly f!75|) . The "suitable boundary theory" is thus exactly the chiral theory 
whose anomaly we obtained in the previous section. 

In contrast to C[T], the Chern-Simons functional C[u] is reparametrization invariant, but 
it fails by boundary terms to be invariant under rotations (or Lorentz transformations) of 
the frame field: 

e a -> = e b O b ai 
u\ -> (co°) a b = {0- l )\uj c d O d b + (0- l ) a c dO c b . (92) 

20 



To obtain the energy momentum flow associated with C[u] we should remember that u is 
linked to the metric through the torsion-free condition 



de* a + oo\ A e* b = 0. (93) 

and through g^ v = r\ ab e*^e* b . We therefore define a tensor T bc and its contravariant version 
rpda _ ^<*^acj^ var yi n g vielbein: 

5S cS = J d n x^g (^j 6e» 

= J d n x^ {T bcV ca e; b ) 5e£. 

= j d n x^T da 5e da . (94) 

The last line introduces the useful quantity. 5ed a = f]db€* b 5e^. As defined, there is no 
immediate reason for T bc to be symmetric. However when the functional S is invariant 
under an infinitesimal local rotation 8e^ = e^9 b a , we have 

= 5S cS 

= J d n x^T bcV ca e* b e^ d a 

d n x^T bcV ca e b a 



/• 
/■ 



d n x^T aa 9 da . 

Since 9 da is an arbitrary skew symmetric matrix, we see that T da = T ad . Accepting this 
symmetry, we can now set 



5S eS = \f drxy/gT^efieZ + TpefSeZ) 
= \J d n x^T aP (e^ ca 5e a a +etv 
= \J d n x^T a p5g a P. 



Here T a p = e b a e c pT bc . Thus, for rotation invariant actions, the vielbein variation leads to the 
same energy- momentum tensor as Hilbert's metric variation. 
Now we have 

SC[co] = —[ tr{5tuR} + — [ tr {5udu}. (95) 

487T J M 967T J dM 
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and we can use 







+(V fc 5eji - 


- Vi5e jk ) 


-(Vi5e k j - 


- VjSe k i)} 



to compute T a &. We do not have to perform this rather tedious computation, however. We 
know that the variations of C[T] and C[u] differ only by boundary terms. The bulk energy- 
momentum tensors for the two actions must therefore coincide. The boundary variations 
will differ though. Because C[u] is reparametrization invariant, the Wess-Zumino term 

W[lo,0} = C[lo°} -C[u] (96) 

that together with C[u>] gives the rotation and reparametrization invariant action C[u; ], 
must give rise to a conserved boundary-theory energy- momentum tensor T^ z . This tensor 
must also be covariant under co-ordinate changes, but will not be symmetric. There is only 
one possibility — the frame field version of (1761) : 

T$ z = f ab = f ab - — — —e ab R. (97) 

The contribution X ab that comes from the boundary part of the variation of C[u] will then 
repair the asymmetry. This contribution is easily computed, and is 

X ab = _± l _ e ab R /gg) 

The net effect is that we get the same boundary-theory energy-momentum tensor T^ u = 
T^ z + X^, and the same anomaly equation, independent of whether we write the gravita- 
tional Chern-Simons function in terms of T or in terms of u. The only difference between 
the two formulations lies in the manner in which the boundary energy-momentum is appor- 
tioned between the bulk Chern-Simons contribution X^ v and the boundary Wess-Zumino 
part T^ z . 



VI. CONCLUSIONS 



We have seen that it is most likely that the thermal Hall currents on the surface of 
topological insulators are confined to one dimensional domain walls, and cannot flow in 
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the two-dimensional surface. To confirm this idea we computed the energy-momentum flux 
associated with a gravitational Chern-Simons term in the boundary effective action. We 
found that the energy-momentum flux is proportional to gradients of the Ricci curvature, 
and therefore needs tidal forces to be non-zero. We related this flux to the gravitational 
anomaly experienced by modes localized on one-dimensional domain walls within the surface, 
and showed that this anomaly takes the same covariant form independently of whether we 
write the gravitational Chern-Simons functional in terms of the Christoffel symbol T or the 
spin connection u. 
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